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We present a X-ray spectroscopic study of the bright Compton-thick Seyfert 2s 
NGC 1068 and Circinus Galaxy, performed with BeppoSAX. Matt et al. (1997 and 
1998) interpreted the spectrum above 4 keV as the superposition of Compton reflection 
and warm plasma scattering of the nuclear radiation. When this continuum is extrap- 
olated downwards to 0.1 keV, further components arise. The NGC 1068 spectrum is 
rich in emission lines, mainly due to Kq, transitions of He-like elements from oxygen 
to iron, plus a fluorescent line from neutral iron. If the ionized lines originate in 
the warm scatterer, its thermal and ionization structure must be complex. From the 
continuum and line properties, we estimate a column density, Nwarm, of the warm scat- 
terer less than a few x 10^^ cm~^. In Circinus Galaxy, the absence of highly ionized iron 
is consistent with a scattering medium with Ux ~ 5 and Nwarm ~ a fewxlO^^ cm~^. 
In both cases the neutral iron line is most naturally explained as fluorescence in the 
medium responsible for the Compton reflection continuum. In NGC 1068 an optically 
thin plasma emission with kT ~500 eV and strongly sub-solar metallicity is required, 
while such a component is only marginal in Circinus Galaxy. We tentatively identify 
this component as emission of diffuse hot gas in the nuclear starbursts. Possible causes 
for the metal depletion are discussed. 

Key words: Line: formation - galaxies: individual: NGC 1068, Circinus Galaxy - 
galaxies: active - X-rays: galaxies 



1 INTRODUCTION 

In the last fifteen years, a wide consensus has gathered 
around the idea that the nucleus and the Broad Line Re- 
gion (BLR) in Seyfert 2s are hidden to us by intervening 
matter along the line of sight (see Antonucci 1993 for a re- 
view on unification models). While column densities « 10^^- 
lO^* cm~^ were measured by Ginga in most Seyfert 2s 
(Awaki et al. 1991; Smith & Done 1996), so confirming this 
hypothesis, no significant absorption in excess of the Galac- 
tic contribution was detected in some objects. A key to ex- 
plain the nature of these apparently odd sources came from 
the evidence that they exhibit intense iron lines, strongly 



suggesting that the observed X-ray emission is due to scat- 
tering of the nuclear radiation, which in turn is completely 
obscured. In fact, if the column density of the absorbing 
matter exceeds 10^"* cm~^), the nucleus is invisible up to 
at least 10 keV (and then no absorption would be directly 
observed). At least two possible candidates exist as refiec- 
tors: the inner side of the molecular torus ("cold refiector") 
and the hot plasma responsible for the scattering of optical 
broad lines ("warm mirror" or "warm reflector"). 

In the former case, the X-ray spectrum should be dom- 
inated by a "bare" Compton-reflection component (George 
& Fabian 1991; GhiseUini et al. 1994; KroHk et al. 1994) and 
an intense (equivalent width EW ^ 1 keV) fluorescent 
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MBF96 hereinafter). This is indeed what has been observed 
in a handful of objects so far: NGC 1068 (the archetypical 
Seyfert 2 galaxy: Ueno et al. 1994; Iwasawa et al. 1997; Matt 
et al. 1997a, M97 hereinafter), NGC 6240 (Iwasawa & Co- 
mastri 1998); Circinus Galaxy (Matt et al. 1996b, Matt et 
al. 1999; M99 hereinafter); NGC 7674 (Malaguti et al. 1998); 
Mkn 3 (Turner et al. 1997a; Cappi et al. 1999). These objects 
constitute the set of so-called "Compton-thick" sources. Re- 
cent results from a BeppoSAX study of an optically-selected 
sample suggest that they are far more common than previ- 
ously thought due to selection biases in the Seyfert 2 samples 
observed in X-rays before (Maiolino et al. 1998). 

The X-ray continuum scattered by the warm material 
is instead a fainter replica of the nuclear non-thermal con- 
tinuum (if self-absorption effects are not important), plus 
a set of emission lines from highly ionized heavy elements, 
which are produced by fluorescence/recombination and/or 
resonant scattering in the ionized matter (MBF96; Netzer 
1996; Netzer et al. 1998, NTG98 hereinafter). 

The X-ray observations, performed so far, suggest that 
both mechanisms are usually at work. Soft excesses above 
the extrapolation of the intermediate X-ray 2-10 keV power- 
laws have been revealed in 21 out of 25 Seyfert 2s observed 
by ASCA (Turner et al. 1997b), and in all "Compton-thick" 
sources (Turner et al. 1997b; Iwasawa & Comastri 1998; 
Matt et al. 1996b). Emission lines from He-, H-like iron and 
lighter elements have been detected in 4 out of 7 Compton- 
thick sources (Matt et al. 1996b; Turner et al. 1997b). Scat- 
tering models can generally account well for these soft ex- 
cesses. However, a description of these features in terms 
of optically thin thermal emission often yields statistically 
comparable results. This is not surprising. Some studies in- 
dicate that the host galaxies of Seyfert 2s have high levels of 
star formation (Maiolino & Rieke 1995) and several exam- 
ples of nuclear starbursts occurring in Seyfert 2s are known. 
The average contribution of starbursts to the 0.5-4.5 keV 
flux was estimated to be 60% in the sample of Seyfert 2s 
observed by ASCA (Turner et al. 1997b), although admit- 
tedly ASCA could not disentangle unambiguously the star- 
burst from the scattered emission. In only one case a signifi- 
cant contribution from starburst could be ruled out (Mkn 3; 
Turner et al. 1997a). 

High spatial resolution is required to separate the con- 
tribution of different spectral components and to perform 
spatially-resolved spectroscopy. Hopefully, this capability 
will be provided in the near future by the detectors on board 
Chandra and XMM. In the meantime, a wide enough spec- 
tral coverage can be exploited to separate the cold and warm 
reflection components above 3 keV and to quantify any con- 
tribution to the soft X-ray emission from a starburst com- 
ponent. The Italian-Dutch satellite BeppoSAX (Boella et al. 
1997a), whose scientiflc payload covers the energy range be- 
tween 0.1 and 200 keV, is the presently operating mission 
best suited for this purpose. For example it has allowed for 
the first time the measurement of the relative contributions 
of the cold and warm mirrors in NGC 1068, confirming the 
complex nature of the reflector, flrst suggested by ASCA 
spectroscopy of the iron line complex (Iwasawa et al. 1997) . 
In this paper we will focus on the soft X-ray properties of 
NGC 1068 and Circinus Galaxy, which were observed by 



thick Seyfert 2s. 

Both sources are very pertinent to the above discussion. 
In both of them, broad optical lines have been detected in 
polarized light (the degree of polarization being ~16% in 
NGC 1068, Antonucci & Miller 1985; and ~2% in Circinus 
Galaxy, Oliva et al. 1998). Both galaxies are site of strong 
starburst activity. In NGC 1068 the bulk of this activity 
is concentrated in a ring of approximately 1 kpc size [i.e.: 
15-16" at the distance of 14.4 Mpc (TuUy 1988)], which pro- 
trudes a bar towards the nucleus (Scoville 1988). The HRI 
image revealed that half of the X-ray emission at 0.8 keV 
comes from an extended region of ~ 13 kpc scale (Wilson et 
al. 1992). It is worth remembering that NGC 1068 hosts a 
"water maser" source (Claussen et al. 1984), probably coin- 
cident with the inner region of a dusty, nearly edge-on torus 
at a distance of ~0.4 pc from the nucleus (Greenhill et al. 
1996). 

Circinus Galaxy also exhibits a strong and variable H2 
maser emission (Greenhill et al. 1997) , which instead may be 
produced in a ^10 AU Keplerian disk around the nucleus. 
Enhanced star formation activity in the shape of a ring is 
present on a ~200 pc scale (corresponding to ~10" at the 
distance of 4 Mpc; Marconi et al. 1994). Although the total 
luminosity of the starburst within a few hundred pc is com- 
parable to the intrinsic luminosity of the Seyfert nucleus, 
only 2% of it is radiated within the inner 12 pc (Maiolino et 
al. 1998a). 



2 OBSERVATIONS, DATA REDUCTION AND 
PREPARATION 

The Itahan-Dutch satellite BeppoSAX (Boella et al. 1997a) 
carries four co-aligned Narrow Field Instruments. Two of 
them are gas scintillation proportional counters with imag- 
ing capabilities: the Low Energy Concentrator Spectrometer 
(LEGS, 0.1-10 keV, Parmar et al. 1997) and the Medium 
Energy Concentrator Spectrometer (MEGS, 1.8-10.5 keV, 
Boella et al. 1997b). They have an energy resolution of 
~8% at 6 keV (MEGS) and of ~4% at 1 keV (LEGS). The 
other two instruments are coUimated detectors, mounted 
on a rocking system to allow a continuous monitoring of 
the background: the High Pressure Gas Scintillator Propor- 
tional Counter (HPGSPC, 4-120 keV, Manzo et al. 1997) 
and the Phoswitch Detector System (PDS, 13-200 keV). 
The HPGSPC is tuned for spectroscopy of bright sources 
with good energy resolution, while the PDS possesses an un- 
precedented sensitivity in its energy bandpass. Only LEGS, 
MEGS and PDS data will be considered in this paper, since 
the HPGSPC failed to detect both sources. 

The NGC 1068 observation was interrupted after ~70% 
of the scheduled exposure time and completed about one 
year later. The log of the observations is reported in Table |l| 
Data were reduced according to the same prescriptions as in 
M97. The only difference regards the PDS data, on which 
a crystal temperature dependent Rise Time threshold was 
applied during the screening procedure. By reducing the in- 
strumental background by ~ 50%, this method allowed an 
improvement in the S/N ratio from 4.5 to 8a". We assume 
hereinafter that this also implies a decrease by a factor 0.92 
of the effective area in comparison to the publicly available 
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lable 1. UcppobAA observations log. icxp and Uit are me total eltective exposure time and count rate in tne U.i-4 keV, i.o-iU.5 keV 
and 13-200 keV bands for the LEGS, MEGS and PDS, respectively. 



Source 


Start time (UTG) 


End time (UTG) 


rpLECS 

cxp 

(s) 


CrLECS 

(B-l) 


rpMECS 

cxp 

(s) 


CRMECS 
(s-i) 


rpPDS 

cxp 

(s) 


GRPDS 


NGC 1068 (1) 
NGC 1068 (2) 
CG 


30/12/96 08:47:25 
11/01/98 09:52:36 
13/3/98 06:33:01 


03/01/97 05:27:50 
12/01/98 08:07:50 
17/3/98 06:21:05 


61497 
15408 
83853 


0.1163 ±0.0014 

0.110 ±0.03 
0.0162 ± 0.0005 


100150 
37331 
137700 


0.0872 ±0.0011 
0.0664 ±0.0015 
0.0962 ± 0.0009 


62493 
17657 
63200 


0.18 ±0.02 
0.15 ±0.06 
2.01 ±0.04 



matrix (September 1997 release, Fiore et al. 1998), which 
were used throughout. Background spectra for the imaging 
instruments were extracted from blank field event files, using 
the same region in detector coordinates where the sources 
lie. PDS background subtracted spectra were obtained by 
plain subtraction of the "off-source" from the "on-source" 
ones. In all spectral fits, multiplicative factors were included 
to account for the absolute flux cross calibration misalign- 
ment between the BeppoSAX detectors (Grandi et al. 1997; 
Haardt et al. 1998). The LEGS vs. MEGS factor was left 
free to vary as a free parameter in the fits, and turned out 
to be ~0.9. The PDS vs. MEGS factor was fixed to 0.75. The 
results are not substantially affected by the residual ^ 5% 
uncertainty on the last parameter (Fiore et al. 1998). 

In the following, errors are at the 90% confidence level 
for one interesting parameter {i.e.: Ax^ ~ 2.71); energies are 
quoted in the source rest frame; Ho — 50 km s~^ Mpc~^, 
and cosmic abundance after Anders & Grevesse (1989) are 
assumed. 



3 NGC 1068 

Source spectra for NGG 1068 were extracted from circular 
regions of 8' and 6' radius for the LEGS and MEGS, respec- 
tively. The spectral analysis was performed in two steps. 
First, we re-analyzed the high-energy (i.e. above 4 keV) 
spectrum, to check the results of M97 after the improved 
PDS data screening algorithm and also to include the 1998 
data set. Once the high-energy continuum shape was es- 
timated, we tackled the task of describing the broadband 
spectrum. We have not found any difference in the spec- 
tral parameters larger than the statistical uncertainties be- 
tween the 1996 and 1998 observations in any of the spectral 
model employed in this paper. The spectral analysis have 
therefore been performed on their unionQ ("total" dataset 
hereinafter) . 

3.1 E > 4 keV 

In the analysis of the NGG 1068 spectra, we have fol- 
lowed the same approach as in M97. The MEGS and PDS 
spectra above 4 keV have been fitted simultaneously with 

* The analysis of the total spectrum has been performed on the 
LEGS and PDS summed spectra, while the MECS spectra of the 
two observations have been maintained separated and fitted si- 
multaneously given the different number of operative MECS units 
(three versus two) between the two observations. 



a model composed by a Gompton reflection continuum 
(model pexrav in Xspec, Magdziarz & Zdziarski 1995) and 
a power- law (describing the warm reflector component), 
both of them absorbed by the Galactic column density 
(Nh = 3.1 X 10^" cm-2, Dickey & Lockman 1990). The pho- 
ton index of the power-law and of the illuminating primary 
continuum, which gives rise to the Gompton reflected com- 
ponent, were tied together. No cut-off has been assumed 
in the primary component. To describe the iron line com- 
plex we first tried a system of three narrow (i.e.: intrinsic 
width a equal to 0) Gaussian lines; the centroid energy E 
of one of these components was held fixed to 6.4 keV (i.e.; 
fluorescence from neutral or mildly ionized iron) while the 
other two were left free to vary. The addition of a further 
narrow line yields a Ax^ = 8.3 for two degrees of freedom, 
significant at more than 97.5% level of confidence, giving 
a line energy of ~ 8.1 keV. The quality of the fit is very 
good (x^ = 165.7/179 dof) and the photon spectral index 
(r ~ 2.0 ± 0.2) is consistent with that typically observed 
in Seyfert 1 galaxies (Nandra & Pounds 1994; Nandra et 
al. 1997). The best-fit solution requires the centroid en- 
ergy of the Ka ionized iron lines to be ~ 6.561q q7 keV 
(corresponding to Fexx+j") and 6.96+o;(')8 keV (Fe xxvi), 
with Equivalent Widths (EW) against the total continuum 
of 740 and 820 eV, respectively. The 8.1 keV line has a cen- 
troid energy consistent with the Kfs transition from H-like 
iron (Ec = 8.1 ± 0.2 keV) and its ratio with the 6.96 keV 
line is 0.18 ± 0.12. The MEGS does not have enough energy 
resolution to allow a totally unambiguous deconvolution of 
the iron line complex in this source. Gomparably good fits 
(albeit formally slightly worse) can be obtained also if one 
assumes a system of Fel+6.61 keV+6.86 keV lines (as sug- 
gested by Iwasawa et al. 1997; = 173.6/180 dof) or of 
Fel+Fexxv+Fexxvi (x^ = 171.5/180 dof). The 8.1 keV line 
is required at comparable statistical significance in all the 
above deconvolutions. The relative weight of the iron 
components varies significantly in the three models. In the 
following we will assume the last one as our baseline as it is 
the most physically reasonable. The properties of the lines 
in this scenario are reported in Table M. 

In Figure ^ the best-fit model and residuals photon 
spectrum is shown for the total dataset. It will be referred 
to as the "double-reflector" model hereinafter, whereas 
the double-reflector continuum will refer to the above 
model without the iron emission line complex. 
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Figure 1. NGC 1068 spectra {upper panel) and residuals in units 
of standard deviations {lower panel) when the best-fit double- 
reflector model is applied to the MECS/PDS spectra above 4 keV 
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Figure 2. Extrapolation of the high-energy best— fit double- 
reflector model {solid line) in the LEGS energy bandpass. Note 
the huge soft excess observed {crosses). The extrapolation of a 
power— law having the same flux at 4 keV but spectral index 
r = 1.2 (as observed by Ueno et al. (1994, dotted line) is also 
shown for comparison 



3.2 Broadband spectrum 

In Figure ^, the extrapolation of the double-reflector model 
in the LEGS band (below 4 keV) clearly shows the presence 
of a huge soft excess. Although this component has already 
been pointed out by several authors in the past (Marshall 
et al. 1993; Ueno et al. 1994), BeppoSAX unprecedented 
broadband coverage allows for the first time a self-consistent 
simultaneous study of both the hard emission and the soft 
excess properties. It is worth stressing that the soft excess 
is present despite the fact that the extrapolated spectrum 
is now much steeper than assumed so far. Previous studies 
have in fact calculated a mean intermediate X-ray spectrum 
and estimated the soft excess above it (cf. Figure g). Now 
we can deconvolve the observed flat spectrum in two physical 
components. 

We have modeled the soft excess by adding to the 




Energy (keV) 

Figure 3. Data/model ratio when the double reflector model plus 
a thermal optically thin plasma is applied to the LEGS {filled 
circles), MEGS {empty squares, 1996 dataset; empty circles, 1998 
dataset) and PDS spectra (not shown) simultaneously. The labels 
indicate the likely identiflcation of the observed lines (marked at 
the expected energies after taking into account the redshift of the 
source). Each data point has a signal-to- noise ratio > 10 



Table 2. Best-flt results obtained for NGG1068 when an extra 
component (indicated in the flrst row) is added to the double- 
reflector model in order to describe the observed soft excess. 
po = power-law; br = thermal bremsstrahlung; mk = optically 
thin plasma 



Parameter 



po 



br 



mk 



I^hard 
Tsoft 

Nsc/Npcxrav (%) 

kT (eV) 
Az (%) 
xVdof 



5.0 ±0.4 

o 21+013 
< 1.2 



3.1 ±0.3 

2 12+01'' 



7±t 



410±t 



513.0/442 453.7/440 



2.4 ±0.4 
2.13 ± 0.17 

7±4 
440 ± 50 
2.8 ± 1.6 
451.3/441 



double-reflector model a thermal plasma emission (model 
mekal in XSPEC), or a bremsstrahlung or a power-law. Al- 
though the quality of the flt improves dramatically in each 
case, local wiggles can still be seen both in the LEGS and 
in the MEGS spectrum (see Figure suggesting that lo- 
calized emission features are required by the data as well. 
We have therefore followed the approach of adding narrow 
Gaussian lines, until statistically required by the flt at the 
99% signiflcance level. This yields a system of flve soft X-ray 
lines in the range 0.5-2.5 keV. The best-flt continuum and 
line parameters are shown in Table ^ and |^, respectively. 

When the soft excess is modeled with an optically thin 
plasma, the is rather good {xt ~ 1.02). The addition 
of a further thermal component with a different tempera- 
ture is not statistically required by the data {Ax^ < 0.1). 
Note that the nuclear power-law index remains consistent 
with that typically observed in Seyfert 1 galaxies, although 
formally slightly steeper than the value obtained when the 
high-energy spectrum alone is analyzed. The thermal com- 
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Energy (keV) Energy (keV) 

Figure 4. Left panel: Spectra and residuals in units of standard deviations when the best-fit model (double reflector model plus optically 
thin plasma) is applied to the broadband NGC 1068 spectrum (only LEGS and MEGS data are shown for clarity). Each data point 
has a signal-to-noise ratio > 10. Right panel: inferred best fit model. The Solid line represents the total spectrum, the dashed lines the 
continuum spectral components and the dotted lines the emission lines 



ponent best-fit abundance is strongly sub-solar (~ 3%), the 
temperature is ~440 eV and the unabsorbed total luminosity 
is -g-SxlO"*^ erg s'^ (0.1-10 keV). The best-fit line centroid 
energies are 0.54, 0.95, 1.32, 1.88, 2.45 keV, and are consis- 
tent (within the statistical uncertainties) with fluores- 
cence of He-like oxygen, neon, magnesium, silicon and sul- 
phur {i.e.: Ovil, Neix, Mgxi, Sixil, Sxv) (see Table |). The 
best-fit results and model are shown in Figure ^ The ther- 
mal component contributes substantially only to the iron-L 
complex. The soft X-ray lines EW against the total contin- 
uum are in the range 70-180 eV, while against the scattered 
continuum are in the range 130-340 eV. 

The total observed fluxes in the 0.1-2 and 2-10 keV 
bands are 1.11 and 0.52 x 10~^^ erg cm~'^ s~^ in the 
0.1-2 and 2-10 keV bands, respectively. They corre- 
spond to unabsorbed rest frame luminosities of 1.31 and 
0.38 X 10*^ erg s"^ The relative fluxes (2-10 keV band) of 
the Compton reflection and scattering components are 1.9 
and 2.2 x 10~^^ erg s~^ cm~'^, corresponding to unabsorbed 
luminosities of 1.34 and 1.55 x lO**^ erg s~^, respectively. 

The abundance of the best-fit mekal model in the above 
scenario is very low and the contribution to the observed soft 
X-ray emission lines is therefore in most cases negligible. 
We, therefore, repeated the fit, assuming this time that the 
soft excess is simply due to a featureless continuum compo- 
nent, either a thermal bremsstrahlung or a power law (even 
though these models could not be physically plausible) . The 
best-fit results and parameters are shown again in Tables ^ 
and ^ In the thermal bremsstrahlung case, the main differ- 
ence is that a line at E ~ 0.78 keV is significantly required, 
to account for the missing mekal contribution to the iron-L 
complex (see Figure^). This new line has an EW of ~100 eV 
(380 eV) against the total (scattered) continuum. Not sur- 
prisingly, the other best fit parameters are very close to those 
obtained in the optically thin plasma scenario. By contrast, 
a much worse fit is obtained if instead a power-law is used 
to model the soft excess (only an upper limit on the scatter- 
ing fraction can be set in this case). The fit is unacceptable 
at > 98% level of confidence {xl = 1.16 for 442 dof) and 



the steeper index obtained for the intrinsic nuclear emission 
(r ~ 2.3) produces a systematic underestimate of the PDS 
counts. 

It is worth noticing that an unacceptable fit is obtained 
(X^ = 900.1/448 dof) if we abandon the double- refiector sce- 
nario and assume that the soft excess (above a bare Comp- 
ton refiection) and the lines from ionized species are ac- 
counted for by several single-temperatures optically thin re- 
gions. Similar results are obtained if one adopts a single 
multi-temperature emission plasma. This provides a further 
indirect confirmation of the validity of the high-energy spec- 
tral deconvolution. 



4 CIRCINUS GALAXY 
4.1 BeppoSAX data 

The ASCA 0.5-10 keV spectrum of the Circinus Galaxy is 
well fitted by a bare Compton-reflection component, a soft 
excess and a system of emission lines (Matt et al. 1996b). 
The BeppoSAX observation has permitted to discover that 
at energies above ~10 keV the primary continuum emerges, 
suggesting that the nuclear region is seen through a screen 
of absorbing matter with Nh ~ 4 x 10^'' cm"'^. The reader 
is referred to M99 for a detailed study of the high-energy 
spectrum. In the following we will focus on the soft excess 
and emission line system. Only LEGS and MEGS data will 
be considered in the following, since PDS data are largely 
affected by the transmitted component, whose contribution 
is instead negligible in the 0.1-10 keV band. For the data re- 
duction see M99. Source spectra have been extracted from 
circular regions of 2' radius around the source image cen- 
troid, to avoid contamination from a serendipitous source in 
the field of view (M99). 

The available statistics below 2 keV is limited by the rel- 
atively large Galactic absorption (Nn.Gai ~ 3 x 10^^ cm~^. 
Dickey & Lockman 1990). There is therefore less room than 
in NGG 1068 for deconvolving the soft X-ray spectral com- 
plexity and/or to perform a detailed line spectroscopy 
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lable a. A-ray lines detected m iNLrU iUbo. i^entroid energies are 
in keV, EW in eV. EW^^^ is the equivalent width against the total 
continuum, EgJ.' against the warm scattered continuum only, 
against the Compton-reflection continuum only. The third column 
indicates the most likely identification. In brackets the line fluxes 
are quoted in units of photons cm~^ s~^. The properties of the 
iron line complex are here given for the Fel+Fexxv+Fexxvi plus 
the ~8.1 keV line deconvolution. All lines refer to Ka transitions, 
unless otherwise specified. 



Parameter 


br 


mk 


Identification 


E(i) 


n 5Q + n 04 


n 54 + 04 


Ovil 




Mil 


80±|0 


70 ±30 


(1.0 X 10- 






340±?|0 


310 ±140 






E(2) 


n 70+0.06 

'^•"'='=0.04 

100 ± 30 


O.ygt 

< 50 


Fe-L 




EwiV 


380 ± 110 


< 180 






E(3) 

Mil 


n q>;+0.04 

150±50 
480±l™ 


yj.jij _i_ yj.yj-j 

90 ±30 
320 ± 100 


Neix 
(3.3 X 10- 




E(4) 


i.oznig 05 
80 ±30 


1 S2 + 05 

-L.Oii _l_ yj.yj'j 

70 ±30 


Mgxi 
(8.7 X 10- 




EW(^' 


180 ± 70 


170 ± 70 






E(5) 

tot 


1.88±0:03 
180±|0 


1.88 ±0.03 
170 ± 40 


Sixiii 
(6.9 X 10" 




EW^c' 


280±™ 


280 ± 60 






E(6) 


2 46+0.03 

Z.^IDICq 04 

160 ± 40 


2.45 ± 0.04 
160 ± 40 


Sxv 
(3.1 X 10- 




EWL'^ 


220 ± 50 


210 ±50 






E(7) 


6.7t 


Fexxv 




EW*:) 


1000 ± 200 


(5.4 X 10- 






3000 ± 600 






E(8) 


6.96t 


Fexxvi 




EW(^) 


530±;« 


(2.5 X 10- 




EWL" 


1500±|00 






E(9) 


8.1 ±0.2 


Fexxvi (K^) 


Mil 


230 ± 130 


(7.0 X 10- 


-6 


EW^?' 


600 ± 300 






E(iO) 


e 


.4t 


Fei-xvi 




EW(^») 


1000 ± 140 


(5.5 X 10" 




EwfJ"' 


1600 ± 200 








SiXIII ArXVII 

sxv 



Fe,K« 



Ni 



1 ^ 



Energy (keV) 

Figure 5. Spectra (upper panel) and residuals in units of stan- 
dard deviations {lower panel), if the double-refiector contin- 
uum + one Gaussian iron fluorescent K„ line model is applied to 
the LEGS {filled circles) and MEGS {crosses) spectra of Circinus 
Galaxy. Most likely identification of the residual line features are 
labeled. All lines refers to Kq, transitions if not otherwise stated 



Table 4. Same as Table 3 for Gircinus Galaxy. 



t fixed 



E(i) 
EW. 

EW; 

E(2) 
EW. 

EW; 

E(3) 
EW. 

EW; 

E(4) 
EW. 
EW, 

E(5) 
EW. 
EW, 

E(6) 
EW. 

EW; 

EW, 



(1) 

tot 

(1) 



(2) 
tot 
(2) 



(3) 
tot 
(3) 



(4) 
tot 
(4) 



(5) 
tot 
(5) 



(6) 
tot 
(6) 



i-87±g:8| 

120±fO 
160±i° 

2 43+0.05 

190 ± 50 
260 ± 70 

3 14+0.11 
•-'•^^='=0.09 

90±|0 
140±60 

lO.012 
^0.013 

2250 ± 100 
2850 ± 130 

7.08 ±0.06 
500±|0 
640±i;o 

180 ± 90 
700 ± 300 
240 ± 120 



1.446±" 



Sixiii 
(3.4 X 10-5) 

Sxv 
(4.0 X 10-5) 

Arxvn 
(1.5 X 10-5) 

Fexvn-xix 
(3.19 X lO-'') 

"Neutral" iron (K^) 
(6.1 X 10-5) 

Fexxv (or Ni K^) 
(1.6 X 10-5) 



A rather bad fit is obtained if a double-reflector 
continuum plus an iron fluorescent line is used 

(X^ = 594.2/78 dof, see Figure H). Most of the residuals are 
due to several unaccounted emission lines. We tiien added 
narrow lines one-by-one, until required in terms of F-test 
statistics at 99% confldence level. Eventually, six lines were 
needed, with best-flt centroid energies of 1.87, 2.44, 3.14, 
6.45, 7.09, and 7.9 keV (see Table |). The most hkely iden- 
tifications for the first five lines are Kq. fluorescence of Sixill, 
Sxv, ArXVII, Fexvii-xix, and Fe<xvii Kp, respectively. 
Most of them had already been measured in the ASCA ob- 



servation of Circinus Galaxy, with comparable EW (Matt et 
al. 1996b). It is worth noticing that the ~7.09 keV line is in- 
consistent, given the uncertainties, with being produced by 
fluorescence of Fexxvi, resolving an ambiguity that was still 
present in the ASCA data. The huge (EW ~ 2.3 keV) hue 
from neutral or mildly ionized iron is well consistent with 
being produced in the same cold reflector, whose emission 
dominates the 2-10 keV spectrum (Matt et al. 1996b). 

The identification of the ~7.9 keV line is not straight- 
forward. In principle, it could be associated with K0 fluo- 
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expect a substantial contribution from the Kq, line of this 
ion as well. Only a 90% upper limit of 52 eV can be set 
on the EW of a 6.7 keV line (against the total continuum), 
implying a K/s/Ka intensity ratio >0.22. Alternatively, the 
line can be due to Kq fluorescence of nickel. The implied ion- 
ization state is still rather high (>Nixx). Again, one should 
expect for consistency a significant line contribution from 
Ka fluorescence of H-like iron at 6.96 keV. The 6.45 keV Ka 
line profile is narrow, the 90% upper limit on its intrinsic 
width being only 60 eV. This excludes in principle any signif- 
icant blending with lines produced by ions more ionized than 
Fexx. The 90% upper limit on the EW (against the scatter- 
ing continuum only) of a 6.96 keV narrow line is 530 eV, 
implying an implausible value of the ionized nickel versus 
ionized iron line ratio ^ 0.65. The accuracy in the MEGS 
instrumental gain reconstruction is ^ 0.3%, and therefore 
comparable with or even slightly worst than the statistical 
uncertainties, but not enough to confuse the observed Kq 
iron line with that expected from an H-like stage. Future 
high resolution, large sensitivity observations are needed to 
allow a more definite conclusion on the identification of such 
a feature. 

The ratio between the neutral/mildly ionized iron Kp 
and Kc, intensities is 0.19 ± 0.04, therefore slightly higher 
than expected from neutral matter (0.11, Kikoin 1976). The 
Kp line energy overlaps with the strong absorption edge as- 
sociated with the Compton reflection component. Any small 
change in the detailed shape of the edge (due, for example, to 
a low degree of ionization of the reflector that might smear 
the edge shape) could have a signiflcant effect on the Kfs 
line intensity. We have tested this possibility using the model 
pexriv in XSPEC (Magdziarz & Zdziarski 1995), which takes 
self-consistently into account the ionization structure of the 
reflector. We have fixed the temperature to the value 10^ K 
to avoid over-fitting the data, thus leaving only the ioniza- 
tion parameter (^) free. The improvement in the quality of 
the fit is very low (Ax^ = 0.8) and the ionization parame- 
ter is very loosely constrained [log(^) = 0.0±i;4]. Nonethe- 
less, the Kp/Ka intensity ratio becomes 0.15±o;o|, consis- 
tent with the expectations. This may be a further piece of ev- 
idence in favor of the idea that the cold reflector in Circinus 
Galaxy is actually mildly ionized. Alternatively, if the iron 
abundance is larger than solar, the intensity of the 7.09 keV 
line could be overestimated to compensate for the missing 
continuum photons above the photoelectric edge of neutral 
iron. Interestingly enough, an iron overabundance by a fac- 
tor of about three is independently suggested by the EW 
of the Kq line (see Sect. 5.1). However, if we leave the iron 
abundance free to vary in the fit, no improvement is ob- 
tained, the K/3 versus Kq intensity ratio remaining basically 
unchanged, and Zpc — 0.9±o;3. 

The ~ 3.14 keV Arxvil line observed by BeppoSAX 
escaped detection in the ASGA data. We stress here that 
the exposure time of the BeppoSAX observation is about 
four times higher than the ASGA one. On the other hand, 
two lines seen by ASGA have not been detected by the LEGS 
(at ~0.8 keV and at ~ 1.33 keV (Mgx-xi)), not surprisingly, 
given the larger ASGA sensitivity at these energies. The cor- 
responding BeppoSAX 90% upper limits on the equivalent 
widths (121 and 53 eV, respectively) are, however, consistent 
with the ASGA values (Matt et al. 1996b). 



iable b. berendipitous sources detected m the KUbAl /titiL ob- 
servation of the sky region containing the Circinus Galaxy 



Source # 


R.A.(J2000) 


Dec.(J2000) 


Count Rate 








(s-i) 


1 


14 13 09.5 


-65 20 14.8 


0.0128 ± 0.0019 


2 


14 13 07.5 


-65 20 06.4 


0.0030 ± 0.0010 


3 


14 13 12.0 


-65 20 07.5 


0.0039 ± 0.0011 


4 


14 12 38.8 


-65 23 30.7 


0.0167 ± 0.0022 



In Figure |6| the spectra, best-fit model (double reflector 
continuum plus emission lines) and residuals are shown in 
Figure]^. As already observed by M99 the best-fit photon in- 
dex of the primary continuum (~ 1.6; see Table ^), is slightly 
fiatter than that typically observed in Seyfert 1 galaxies, al- 
beit not inconsistent given the relatively large uncertainties. 
The total observed fluxes in the 0.1-2 and 2-10 keV en- 
ergy bands are 8.7 x 10~^^ and 1.37 x 10"" erg cm"^ s"\ 
corresponding to unabsorbed luminosities of 2.9 x 10*" and 
1.34 X 10''^ erg s~^, respectively. The 2-10 keV band fluxes 
of the Gompton reflection and scattering components are 6.3 
and 3.0 x 10~^^ erg s~^ cm~^, corresponding to unabsorbed 
luminosities of 6.0 and 2.9 x 10''° erg s~^. 

If a single temperature mekal model is added to the 
best-fit double refiector continuum plus lines model, a 
marginal improvement in the quality of the fit is obtained 
(Ax^ ~ 3.8 for two more free parameters). The Gompton- 
refiection, scattering and line properties are not substan- 
tially affected by the inclusion of this further component. 
The unabsorbed 0.1-10 keV luminosity corresponding to the 
best-fit nominal temperature (~500 eV) is 1.4 x lO**" erg s~^, 
while the abundance is basically unconstrained. 

No significant improvement in the quality of the fit is 
obtained if the optically thin plasma is substituted by a 
bremsstrahlung or if the soft X-ray power-law spectral in- 
dex is left free to vary. Again, neglecting the double-reflector 
scenario and assuming that the whole soft excess and the 
ionized lines are due to a single or multi-temperature opti- 
cally thin plasma yields an unacceptable (98.0/69 dof). 

4.2 HRI data 

Gontrary to NGG 1068, little can be found in the literature 
about the X-ray imaging of the Gircinus Galaxy. For this 
reason, we present here the results of the only ROSAT/HRI 
observation (spatial resolution ~10" Full Width Half Max- 
imum and 2-band spectral resolution in the nominal 0.1- 
2.4 keV energy bandpass), whose data are publicly avail- 
able. The observation was performed from September 14 to 
15, 1995 for a total exposure time of about 4200 seconds. 
Data were retrieved as event list from the HEASARG public 
archive, images were extracted with Xselect and analyzed 
using the Ximage package (Giommi et al. 1991). PHA chan- 
nels between 2 and 12 (approximately 0.1 and 2 keV) were 
used to optimize the signal to noise ratio; only data with a 
good attitude reconstruction were used. 

Four X-ray sources are detected at more than 3a level in 
the HRI field of view (see Table § and also Figure M) . Source 
#1 is very close to the optical position of Gircinus galaxy 
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lable 5. oest-ht results ot the Uircmus galaxy U.l-lU keV hts tor the continuum pajrameters in: (hrst rowj the double reflector continuum 
plus emission lines model and (second row) the same model plus a thermal plasma. 



Model 
mm 


Nh 
10"^" cm^^ 


r 


Nsc /Npexrav 

(%>) 


kT 
fkeVl 


z 




po 

po+mekal 


25 ±7 
80±i8 


1.6 ±0.2 

-'-"='=0.3 


3.8±t^ 


o.5±8:i 


loot 


58.7/67 
54.9/65 



t fixed 




Energy (keV) Energy (keV) 

Figure 6. Left panel: Spectra and residuals in units of standard deviations when the best-fit model (double reflector continuum plus 
emission lines) is applied to the broadband Circinus Galaxy spectrum.Each data point has a signal-to-noise ratio > 10. Right panel: 
inferred best fit model. The Solid line represents the total spectrum, the dashed lines the continuum spectral components and the dotted 
lines the emission lines. 
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Figure 7. Zoom of the central 10 arcmin region of ROSAT/HRI 
field of view containing the Circinus Galaxy field (con- 
tours) , superimposed to an optical plate 



nucleus, allowing the identification of this X-ray source with 
the AGN. Sources #2 and #3 arc 30-40 arcscc away from 
the nuclear region, while source #4 is located about 5 arcmin 
from the galaxy in the S-W direction. This last source was 
also detected by ASCA (Matt et al. 1996b) and BeppoSAX 
(M99). 



No extended emission is detected around the nucleus 
If we subtract the two sources closer to the nucleus, the 
residual nuclear emission is consistent with the PSF of the 
instrument (Davis et al. 1996). We stress, however, that the 
available exposure time is rather low and more statistics is 
needed before reaching a firm conclusion on this point. The 
flux of source #1 is ~ 4.2 x 10~^^ erg cm~^ s~'^, which is 
about half of that observed by BeppoSAX. Sources #2 and 
#3 lie well within the extraction radius used for the Bep- 
poSAX observation and could in principle contaminate the 
soft X-ray flux. If they are not variable, such a contamina- 
tion could be as high as 25%, and therefore account for a 
substantial fraction of the galaxy soft excess. Although the 
possibility that the soft X-ray flux in the BeppoSAX data 
is entirely to one of these serendipitous sources cannot be 
ruled out a prion, it would however require that at least 
one of these sources varies by almost an order of magni- 
tude up to a luminosity level ^1 lO''" erg s~^, while keeping 
coincidentally a F ~ 1.5. Although we cannot reject this hy- 
pothesis on the basis of BeppoSAX (or ASCA) data alone, 
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following. 



5 DISCUSSION 
5.1 X-ray lines 

The BeppoSAX observations of NGC 1068 and Circinus 
galaxy have shown line-rich X-ray spectra. The lines are gen- 
erally associated with He- or H-like ions of elements heavier 
than oxygen, the iron-L complex, and Kq, (and, possibly K/j) 
fluorescent lines from neutral or mildly ionized iron. Simi- 
lar lines had already been observed by BBXRT and ASCA, 
and most often associated with reprocessing of the primary 
radiation in the nuclear environment (Marshall et al. 1993; 
Turner et al. 1996b; Netzer & Turner 1997, M97). The en- 
ergy resolution of the imaging instruments on-board Bep- 
poSAX is not as good as that of the CCD on board ASCA 
and this limits its capability of deconvolving the K-shell 
iron line complex. However, at least two new results emerge 
from the BeppoSAX spectroscopy: i) the detection of emis- 
sion lines around 8 keV, where BeppoSAX instruments can 
take advantage of a better effective area and comparable or 
lower background than the ASCA SIS; ii) the detection of 
a ~0.57 keV line in NGC 1068, which is most likely due to 
the Ka fluorescence from He- like oxygen. The lack of the lat- 
ter feature in previous spectra had forced to assume ad hoc 
oxygen-poor plasma (Marshall et al. 1993; Netzer & Turner 
1997) in modeling the NGC 1068 warm mirror. This problem 
has now been overcome, given the accurate calibration of the 
LEGS instrument on-board BeppoSAX around 0.5 keV (cf. 
Haardt et al. 1998; Orr et al. 1998). 

A common feature of the spectra presented here is that 
the line contribution from any thermal plasma is negligi- 
ble, with the possible exception of the iron-L complex (cf. 
Figure ^ and ^). This has implications on the nature of 
this spectral component which we tentatively associate with 
the nuclear starburst and discuss in the next section. In 
this section, instead, we will focus on the properties of all 
the other lines, which we assume to originate as fluores- 
cence/recombination or resonant scattering from the cold 
and warm reflectors. 

Kct (and K^) fluorescent lines from neutral or mildly 
ionized iron are likely to originate in the cold reflector. 
In Circinus Galaxy, the lack of ionized iron lines allows a 
relatively unambiguous identification and deconvolution of 
the iron emission line complex. The observed Ka EW of 
~2.3 keV is in good agreement with the one measured by 
ASCA (M97) and corresponds to a value of ~2.8 keV if 
measured against the Compton-reflection continuum alone. 
This value, coupled with the relatively low inclination angle 
(i ^ 40°) derived from the hard X-ray analysis (M99), im- 
plies an iron overabundance by at least a factor of 3 (MBF96, 
Matt et al. 1997b). A line is detected as well, with a ratio 
to the Ka line consistent with that expected on theoretical 
grounds, especially when a possible mild ionization of the 
reflecting medium is taken into account. The situation of 
NGC 1068 is far more complex. The limited MEGS energy 
resolution does not allow to unambiguously deconvolve the 
line complex. We have assumed in the above analysis the 
reasonable hypothesis that the line is given by the superposi- 
tion of a "neutral" , a He-like and a H-like component. In this 



reflection continuum is ~ 1.6 keV. Following MBF96, this 
would suggest a high inclination (in agreement with the idea 
that the line comes from the same region as the water maser 
at J ^ 82°), but the large uncertainties on this value prevent 
any conclusion on this point. 

The other lines observed correspond to He-like Ka 
transitions from O, Ne, Mg, Si, S, Fe and H-like Fe. Re- 
cent works have focused on producing self-consistent mod- 
els for the lines expected from a plasma, photoionized by 
a strongly absorbed Seyfert-like continuum. In the follow- 
ing we will mainly refer to the works by Netzer & Turner 
(1997) and NTG98 and deflne the ionization parameter 
Ux as the dimensionless ratio between the ionizing pho- 
ton flux (integrated in the 0.1-10 keV energy band) and 
the electron density (George et al. 1998). The calculations 
used in these papers take into account both the fluores- 
cence/recombination and the resonance scattering lines, the 
latter process being dominant at densities 10^^ cm~^. We 
first consider the more complex case of NGC 1068. The wide 
range of ionized species observed can be hardly reconciled 
with a single-zone, homogeneous plasma. Assuming a den- 
sity NH,warm = 2.5 X 10'^'^ cm~^, all the lines from interme- 
diate elements require 0.5 < Ux < 5, except for the oxygen 
one (cf. Fig. 2 of NTG98). On the other hand, the highly 
ionized iron lines require Ux > 5. At this high level of ion- 
ization, oxygen is likely to be at least in H-like state, if not 
fully stripped. An intense Ovil line is instead still consistent 
with neutral iron, requiring Ux < 0.5. We therefore conclude 
that the warm mirror must be complex and structured. At 
least three components are needed in NGC 1068, unless the 
oxygen line originates in the same cold refiector as the neu- 
tral iron lines. In Circinus Galaxy, the lack of detection of 
oxygen and highly ionized iron lines overcomes these diffi- 
culties, and the observed set of ionized lines can be easily 
explained with a single-zone scatterer with intermediate ion- 
ization parameter. However, while any contribution from a 
strongly ionized scatterer is ruled out, oxygen (and magne- 
sium) lines are likely to have missed detection due to the 
higher Galactic neutral absorption along the line of sight. 
Actually, ASCA revealed a Mgix line, with an EW of about 
90 eV (Matt et al. 1996b), which is consistent with the upper 
limit derived from the BeppoSAX data. 

5.2 Some numbers on the reflectors 

For a geometry such as that adopted by Ghisellini et al. 
(1994), the amount of cold reflection depends mainly on the 
inclination angle, while that of the ionized scattering de- 
pends almost entirely on the optical depth (or column den- 
sity) of the scattering material (MBF96) . From the 2-10 keV 
flux ratio of the warm vs cold reflection, we may therefore 
derive a relation between i and NH,warm, which is shown in 
Fig. I 

In the NGC 1068 case, the source is likely to be ob- 
served at very high inclination angles, as also implied by 
water maser measurements (Greenhill et al. 1996) and by the 
amount of cold reflection continuum measured (M97) . Then 
NH.warm Is coustraiued to be lower than, say, a fewxlO^^ 
cm~'^. An independent estimate of Nn.warm can be derived 
from the equivalent widths of the ionized lines (with respect 
to the warm scattering continuum). These values (~3 keV 
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Figure 8. cost vs. NH^warm relation in NGC 1068 (crosses) and 
Circinus Galaxy {filled circles) 

and ~1.5 keV for the He- and H- like ions, respectively) may 

be compared with Figs. 4, 5 and 6 of MBF96, where the EW 
of these two ions are shown as a function of NH,warm (the 
total line EW decreases with Nn.warm because the resonant 
scattering line becomes more and more optically thick). The 
values in those figures are for solar abundances and fraction 
of ions equal to 1. The value of the EW of the cold iron line 
(with respect to the Compton reflection) implies that the 
iron abundance is roughly solar (MBF96). Assuming typi- 
cal ion fractions of ^ 0.5, the comparison of observed and 
expected EW yields NH.warm values lower than 10^^ cm~^, 
even allowing for turbulent motions in the matter (which 
reduce the line opacity). The EW of the other ionized lines 
are much lower than that of iron, but they are probably cal- 
culated against the wrong continuum and so not much in- 
formation can be gathered from them. In any case, it must 
be stressed that the above calculation is not self consistent 
due to the presence of mildly ionized material responsible 
for the intermediate Z lines; instruments with much larger 
sensitivity and far better energy resolution are necessary to 
make any further progress in understanding the geometry 
and physics of the nuclear cnviroumcnt of NGC 1068. 

The same calculation can be applied to Circinus Gala^xy. 
If we assume that the inclination angle is ^ 40° (M99), the 
density of the warm scatterer has to be higher than a few 
10^^ cm^^. This is in good agreement both with the value 
derived from the Sxv line EW (NTG98) and with the lack 
of detection of highly ionized iron lines (but the 90% upper 
limits on the EW of the He- or H-like iron lines are ~1.5 keV, 
i.e. consistent with practically any column density). The in- 
trinsic luminosity is given by Lx{2 — lOkeV) ~ tA^Jatt ' 
where AO is the solid angle subtended by the visible part 
of the illuminated ionized matter. From the hard X ray 
spectrum, M99 derives Lx{2 - WkeV) ~ lO"^ erg s"\ 
which implies tAQ/Att ~ 0.008. In order to have, say, 
Nh ,warm — 3 X 10^^ cm-^ An/47r should be -0.3, corre- 
sponding to an half-opening angle of the torus of ~65° . 

A NH.warm ~ 10^^ cm~^ scattcrcr should also imprint 



tion of photoionization absorption edges docs not improve 
the quality of the fit in either NGC 1068 or Circinus Galaxy. 
The upper limits on the optical depths of Ovii and Oviii 
photoionization edges are 0.35 and 0.25 in NGC 1068, re- 
spectively. The same limits are 2.7 and 1.5 in the low Galac- 
tic latitude (and therefore highly absorbed) Circinus Galaxy. 
Assuming that the scattering plasma is in the typical condi- 
tions of the "warm absorbers" observed in Seyfert 1 galaxies 
(Reynolds 1997; George et al. 1998), these values correspond 
to hydrogen column densities of NH.warm ^ 5 x lO^'^ cm~^ 
and NH.warm 4 X 10^^ cm~^, broadly in agreement with 
the above diagnostics. 

5.3 On the nature of the soft X-ray continua 

In both Circinus Galaxy and NGC 1068, the soft X-rays 
continuum is most convincingly explained as the superpo- 
sition of two components (even if for Circinus the evidence 
for a second component is marginal) : a scattered power-law, 
which mirrors a fraction fs of the primary nuclear emission; 
and an optically thin plasma emission, which we tentatively 
associate with a nuclear starburst. Indeed, both galaxies ex- 
hibits strong star-forming activity both in optical and IR 
(Scoville 1988; Maiolino et al. 1998). 

Typical temperatures arc ~0.4-0.5 kcV. They are 
slightly lower than those (0.6-0.9 keV) observed in starburst 
galaxies by ASCA (Ptak et al. 1999) and BeppoSAX (Per- 
sic et al. 1998; Cappi 1999). David et al. (1992) studied the 
correlation between far-infrared (FIR) and X-ray luminos- 
ity in starburst galaxies and derived a logarithmic scale law. 
The FIR luminosity of NGC 1068 (calculated assuming the 
formula (1) in David et al. 1992) is 2.2 x lO'"' (Soifer et 
al. 1989); it translates into an expected 0.5-4.5 keV lumi- 
nosity of ~10*^ erg s~^, which is roughly consistent with 
that of the thermal component alone (3.7 x lO**^ erg s^^, 
given also the admittedly wide spread in the correlation. 
For Circinus Galaxy the integrated IR luminosity mear 
sured (4 x lO''^ erg s~^, Moorwood et al. 1996; Maiohno 
et al. 1998) corresponds to an expected X-ray luminosity of 
~ 2.1 X lO**" erg s^^, again broadly consistent with the ob- 
servations. This evidence strengthens the case for the iden- 
tification of this soft thermal component with the starburst. 
However, in Circinus Galaxy we cannot rule out the possi- 
bility that the excess emission is partly or entirely due to 
contaminating X-ray sources in the bulge of the galaxy. 

In NGC 1068 the thermal component abundance is 
more than one order of magnitude lower than solar and in 
agreement with the value derived from the ASCA data (Net- 
zer & Turner 1997). Similar values have been measured in 
starburst galaxies observed by ASCA (Ptak et al. 1998) and 
attributed to the relative weakness of the iron L complex 
lines. When the X-ray spectra of these galaxies are fitted 
with thermal plasmas, and the abundances of iron and a- 
clements are decoupled, the a-clements abundance recovers 
solar or even over-solar values (Persic et al. 1998), whereas 
the iron abundance remains robustly sub-solax (Cappi et al. 
1999). This low metallicity is possibly due to iron depletion 
in warm interstellar clouds surrounded by a dust-rich galac- 
tic environment (Ptak et al. 1997). In principle, dilution by 
featureless continua, which may be provided by an unre- 
solved population of accretion-driven sources, may explain 
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clear starburst environment or in the host galaxy, spatially 
confused with the soft X-ray emission. There is clear evi- 
dence that the nuclear stellar population in NGC 1068 is rel- 
atively young (~5-35 Myr; Davies et al 1998) and, therefore, 
a large contribution by High Mass X-ray Binaries (HMXRB) 
is expected. HMXRB are, however, generally flatter in the 
intermediate X-rays than NGC 1068 and, if they contribute 
significantly to the soft X-ray flux, their summed emission 
should dominate above 1 keV: this is not observed in the 
data. Alternatively, wind accreting Low Mass X-ray Binaries 
(LMXRB) are harder and soft X-ray brighter. In "normal" 
spiral galaxies like our own and M 31 (Trincheri & Fabbiano 
1991), they indeed dominate the total X-ray output, with an 
integrated X-ray luminosity in the range 10^^-10'*^ erg s~^; 
the brightest sources generally cluster around the galactic 
bulge. Of course, the soft X-ray spectrum might well be more 
complex than our simple, statistical-limited two-components 
deconvolution. BeppoSAX and ASCA results, although ob- 
tained with detectors having relatively modest energy res- 
olution, are suggestive of a very complex situation, which 
will be clarified by the forthcoming major X-ray missions 
like Chandra, XMM and Astro-E, whose high resolution soft 
X-ray detectors will surely shed a new light on this subject. 
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